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ABSTRACT 

We present a measurement of the energy carried and dissipated by Alfven 
waves in a polar coronal hole. Alfven waves have been proposed as the energy 
source that heats the corona and drives the solar wind. Previous work has shown 
that line widths decrease with height in coronal holes, which is a signature of 
wave damping, but have been unable to quantify the energy lost by the waves. 
This is because line widths depend on both the non-thermal velocity v^t and the 
ion temperature Ti. We have implemented a means to separate the 7] and v-^t 
contributions using the observation that, at low heights, the waves arc undamped 
and the ion temperatures do not change with height. This enables us to determine 
the amount of energy carried by the waves at low heights, which is proportional to 
Vnt- We find the initial energy flux density present was 6.7±0.7 x 10^ erg cm~^ s~^, 
which is sufficient to heat the coronal hole and acccelerate the solar wind during 
the 2007 - 2009 solar minimum. Additionally, we find that about 85% of this 
energy is dissipated below 1.5 Rq, sufficiently low that thermal conduction can 
transport the energy throughout the coronal hole, heating it and driving the fast 
solar wind. The remaining energy is roughly consistent with what models show 
is needed to provide the extended heating above the sonic point for the fast solar 
wind. We have also studied Ti, which we found to be in the range of 1 - 2 MK, 
depending on the ion species. 

1. Introduction 

One of the major models to describe the heating of the solar corona and the acceleration 
of the solar wind relies on waves to carry the energy. Such wave-driven models have been 
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supp o rted by observations that waves, particularly Alfenic waves (jVan Doorsselaere et al. 



2008 



Goossens et al 



chromosphere ( 



( jTomczyk et al. 



ess et al 



2009[). are observed throughout 



2009 



De Pontieu et al 



20071 ). and into the solar wind ( iBelcher &: Davis 



;he solar atmosphere from the 



2007 



Mcintosh et a. 



19711 ). 



201ll ). to the corona 



However, one difficulty for simple wave-driven models has been that Alfven waves are 
licted to dissipate via viscosil 
the Sun, at about 2-5 Rq (e.g.. 



predicted to dissipate via viscosit y, thermal conductivity, and resistivity relatively far from 



Parker 



1991 



Cranmer 



20021). In order for waves to heat 



the corona they must be damped at much lower heights where heat conduction is more 
efficient. For this reason theories have been developed for how the waves may dissipate more 



quickly. These theories rely on the inhomogeneity of the corona and show, for examp 



the waves can be more strongly damped through phase mixing (IHevvaerts &: Pries 



turbulent cascade (iMatthaeus et al. 



19991 ). or resonant absorption (iGoossens et al 



e, th at 



19831), 



20n|) 



Only recently has incontrov ertible obser v ation al evi dence for dissipatio n of A lfven 



waves at low heights been found. 



Hahn et al. 



torn and lBemporad fc Abbol tom studied 



coronal hole observations and demonstrated that Alfven waves are damped at relatively low 
heights in the corona. In these studies the Alfven waves were observed spectroscopically 
through the non-thermal broadening of optically thin spectral lines. The magnitude 
of the non-thermal broadening is predi c ted to be proportional to the way e amplitdude 



(IBanerjee et al 



1998 



Doyle et al. 



1998 



Moran 



2001 



Banerjee et al. 



20091). For energy to 



be conserv ed the wave ampli tude must increase with 



decreases. 



Hahn et al. 



torn and 



l eight above the Sun as the density 



Bemporad fc Abbo 



(I2OI2I ) found that the line widths 



decrease above about 1.2 Rq. They ruled out systematic errors as the cause of the decrease. 

lat the lin e widths decrease a t larger heights 



This confirmed some 


(Ban 


leriee et al. 


1998 


2008 


)• 



Doyle et al 



1999 



Moran 



2003 



O'Shea et al 



2005 



Dolla fc Solomon 
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In order to determine if waves are indeed responsible for heating the corona and driving 
the solar wind one must quantify both the energy initially present in the waves as well as 
that dissipated by the waves. This has been difficult to determine because the measurement 
of a line width includes contributions both from thermal broadening, which is proportional 
to the ion temperature T;, and from non-thermal broadening, which is proportional to the 
non-thermal velocity v^t- A conventiona l meth od to estimate fnt is to assume some value 



for Tj. For example. 



Bemporad &: Abbd ( |2012| ) assumed the ion temperature was equal to 



the ionization equilibrium formation temperature of the ion emitting the line. However, 
this is not necessarily correct since so me studies have shown that T; ma.y be much greater 



Tu et al. 



1998 



Landi fc Cranmer 



Hahn et al. 



2009 



Hahn et al. 



(I2OI2I ) assumed that Tj = Te 



than the formation temperature (e.g. 
20101). Because Ti is expected to be larger than Te, 
and thereby estimated an upper bound for Vnt- 

Here we present in Section [3] a method to separately determine v^t and T;. The Wnt 
data are presented in Section HI The non-thermal velocity is proportional to the wave 
amplitudes and from it we determine the initial wave energy, the change in the wave energy 
flux density, and the length and time scales over which the waves are damped. These data 
indicate that waves are sufficient to heat the coronal hole and drive the fast solar wind. 
They also provide quantitative constraints for theoretical models of wave damping. In 
Section |5] we present our measurements of ion temperatures and compare them to some 
earlier measurements that found only lower and upper bounds for the temperature. These 
temperature data can be used to test various models for resonant ion heating in the corona. 



2. Observation 



Our data come from four observati ons made wi th the Extreme ult raviolet Imaging 
Spectrometer (EIS: ICulhane et al.l 120071 ) on Hinode ( jKosugi et al.l 120071 ). The observations 
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were made on 2009 April 23 at 12:42, 13:16, 13:50, and 15:17. Each observation was 
30 minutes in duration. For these data the 2" sht was pointed at a polar coronal hole 
at positions relative to the central meridian of X = —14.5", 15.5", 45.4", and 105.6", 



it extendec 



Hahn et al 



rom about 0.95 -Rm to about 



( 120 12l ). but excluding their 



respectively. The height range covered by the si 
1.45 Rq. These data are the same as used by 
observation centered at X = —44.5". We omitted that particular observation as it had a 
larger density scale height at low heights compared to the other pointings, possibly due to 
intervening quiet Sun material. Our results for wave damping, though, are consistent with 
the previous results that included the additional observation. 

The four pointings were averaged together in order to improve the statistical accuracy. 
This was done by first using the standard EIS processing routines to clean the data of 
spikes, warm pixels, and dark current, and calibrate the data. Drif t s in t he wavelength scale 
were then corrected using the method described by 
the data to the same wavelength scale, pixels at the same radius from each of the four 
observations were averaged to create the dataset analyzed. Finally, these data were further 
binned in the vertical direction. For the the analysis described below, where we perform 
a fit to the data at low heights, we have used a binning of 8 pixels per bin (~ 0.01 Rq). 
To extend these results to larger heights, where the intensities are correspondingly much 
smaller, we have used a 32 pixel binning (~ 0.03 Rq). 



Kamio et al. 



fl2010[ ). After aligning 



3. Analysis Method 

3.1. Line Widths 

We fit Gaussian functions to the spectrum in order to derive the line widths. In 
particular, each line was fit with a double Gaussian in order to account for instrument 
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scattered light. The hnes used for various aspe cts of the analysis 



Hahn et al. 



are given in Table [T] 



( I2OI2I ). so here we only briefly 



The fitting procedures are described in detail in 
review the method. 

Stray light is the unshifted spectrum of solar disk emission that is scattered by the 
instrumental optics onto the off-disk data. Because line widths tend to be narrower on the 
disk, scattered light can be a significant source of systematic error at large heights when 
the fraction of real emission is small. To correct for the scattered light we first measured 
line profiles from the portions of our observation that looked at the solar disk. Then we 
constructed a predicted scattered light line profile for each line. For these parameters we 
used the measured line width and centroid position and 2% of the on-disk in tensity. This 



last value is base d on estimates for the magnitude of the scattered light EIS ( lUgarte 



Hahn et al. 



2010 



20121 1. Then, for each position in the off-disk data, a fit was performed using 
a Gaussian with free parameters added to this artificial scattered light profile. This is 
equivalent to subtracting the scattered light profile from the spectrum. 

The measured full width at half maximum AAfwhm of an optically thin spectral line 
depends on ins trumental broadenin g AAinst, the ion temperature T[, and the non-thermal 



velocity v^t as (jPhillips et al. 



20081 ) 



AA 



FWHM 



AALt + 41n(2) 



2fcB^i , 2 



1/2 



Here A is the wavelength of the line, c is the speed of light, is the Boltzmann constant. 



and M is the mass of the ion. We have subtracted the instrumental 



tabulated AAinst values as a function of position along the slit ( lYoung 
width can then be expressed as an effective velocity. 



width using the 



2OI1I). The observed 



VeS 



(2) 



where Vth = \/'ik-QT-J M . This v^e depends on both Ti and 
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3.2. Separating Thermal and Non-thermal Broadening 



DoUa &: SolomonI ( 120081 ) pointed out that the thermal and non-thermal contributions 
can be inferred if two assumptions are made. One can then calculate Vnt(-Ro) at a radius Rq 
using data from another height Ri. The first assumption is that Wth is constant with height 
for each ion emitting the line being studied. This implies that 



^cfr(^i) - ^cfr(^o) = - vURo). 



(3) 



The other assumption is that the waves are undamped. iMoraru (120011 ) has shown that this 
conservation of energy then implies that f^t oc Ue^^^. Since the waves are assumed to be 
undamped we have 



fnt(-Rl 



Vnt{Ro 



n^{Ro) 



1-1/4 



(4) 



Putting it all together, one finds 

Wnt(^o) 



1/2 



no(-Ro) 



-1/2 



(5) 



Doha fc SolomonI ( 120081 ) used this method to determine v^tiRo) by taking an average over 
the results for a fixed i?o while varying Ri. However, they did not observe damping and 
they inferred a quite small v^t, possibly for the reasons we discuss below in Section HI 



Doha fc SolomonI (120081 ) 



The method we use relies on the same assumptions as the 
method, but the application is somewhat different. Here, we use a least squares fit. The 
reason for doing this is that uncertainties in the data can cause large variations in the 
^^nt(-Ro) determined using equation ([5]). A least squares fit takes these uncertainties into 
account implicitly and is more robust to noise. Combining equations ([2|), ([3|), and (jl]), the 
function used in the fit is 



VesiR) = + i4(^o) 



ne{R) 
ne{Ro) 



-1/2 



(6) 



- 8 - 



Here, Vth and fnt(-Ro) are the only free parameters to be determined. As we discuss later, 
we assume that Vth for a given ion is the same for every height, though it can be different 
for each ion. Solving equations (E]) or (E]) requires the ratio ne(-R)/ne(-Ro)- We describe 
below how this is determined. 

The assumption of constant ion temperature is reasonable for low heights. Ions in 
coronal holes are known to be heated, although the precise mechanism has not been 
determ ined. Sorne pos sibilities include ion cyclotron resonance heating by high frequency 



waves fiCranmei 



ion orbits ( IChandran 



2002) and stochastic heating by turbulent fluctuations that disturb the 



LO) 



20101 ). In either case the heating rate is predicted to depend on the 



charge to mass ratio, with minor ions heated more strongly than protons. For this reason, 
at low heights Coulomb collisions with the proto ns are expected to cool the ions and 



maintain them at a relatively steady temperature (iLandi &: Cranmer 



20091 ). Measurements 



estimating upper and lower bounds for Ti at heights of i? < 1.15 Rq have shown that 



is consistent with 



jemg 



(ILandi fc Cranmer 



2009 



constant over this height range, albeit with large uncertainties 



HahnetaL 



2010). Note that each ion may have a different Ti or 



fth, both of which we assume do not change with height. 

It is also reasonable to assume that waves are undamped at low heights, and conse- 
quently fnt oc ^^'^ . This theoretical relation is valid for outw ard propagating waves when 



the s olar wind velocity is much smaller than the Alfven speed (jCranmer fc van Ballegooijen 



20051 ). a condition expected to be met at low heights. Numerous studies have o 



2012 



Bemporad fc Abbo 



met at low 


leights 


(Dovle et al. 


1998 


1 - 



Banerjee et al. 



1998 



Dserved the 



2009; 



Hahn et al. 



2OI2I ). In estimating fnt, these studies have assumed T; to be 



either the ion formation temperature or the electron temperature, but have found the same 
He trend. Thus, this trend is not very sensitive to uncertainties in the magnitude of Ti. 

Based on the above, in the range 1.02 - 1.12 Rq the ion temperatures should be 
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reasonably constant with height; although, Ti may still differ depending on the ion species. 
Also, the upper height of 1.12 Rq is below the point where the waves appear to be damped. 
Thus, it is reasonable to perform the fits to equation ^ over these heights. 

Our analysis also requires an independent measurement of ng. This was obtained 
from the in tensity ratio of the Fe ix 188.50 A and 189.94 A lines using atomic data from 



CHIANTI flDere et al. 



1997 



Landi et al. 



20121 ) . Figure [T] shows the inf erred densities, 



which are typical of densities found in other coronal hole observations ( iWilhelm et al. 



20111 ). At these low heights the solar wind velocity is small and the corona is close to 
hydr ostatic equilibrium. We therefore fit the den sity over the range 1.02 - 1.12 Rq using 



le.g- 



Guhathakurta et al. 



1992 



Doyle et al. 



19991 ) 



ne(-R) = ne(-Ro) exp 



-{R — Rq 



(7) 



HRqR 

where H is the density scale height, which was found to be if = 0.0657 ± 0.0052 Rq. (Here 
and throughout we give the la uncertainties.) The fit is illustrated in figured] Using this 
expression for the density, the ratio ne{R)/ne{Ro) in equation [6] can be rewritten so that 



VesiR) = \ v^^ + vlt{Ro) 



exp 



-{R — Rq 

RRnH 



.-1/2 



This function is useful for our analysis as some of the statistical fluctuations in the 
magnitude of are smoothed out, while retaining the essential description of the density 
variation with height. 



3.3. Uncertainties 

Uncertainties in the line widths v^s and H were propagated into the fitted parameters 
fth and fnt using a Monte Carlo analysis. The observed data for each fit are the single 
inferred value of H and a set of values for Ves as a function of height for a given ion. 
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Normally distributed random variations were added to these input data. The magnitude 
of these variations was set so that the standard deviation of the random numbers added 
was equal to the la uncertainty for each data point used in the fit of equation ([H]) over the 
height range 1.02 - 1.12 Rq (e.g., Figure |2]). The fits were then performed for 1000 different 
variations, with each iteration producing different Vth and fnt. Next we took the mean of 
each parameter and estimated the la uncertainty to be the standard deviation. 

We found that for some of the iterations, the value of i^th for a given ion would imply 
Tj < Te. Although this is clearly a possible fit to the data, the result does not seem 
physically reasonable. At very low heights, below the range of our fit, the density is high 
enough that electrons, protons, and ions should be in equilibrium so that Tj = Tg. At large 
heights the ions are observed to be heated and T; > Te and collisions can be neglected. In 
the range where we perform the fits, the situation lies in-between these two extremes, that 
is the ions are expected to be heated in sonae way, but also to be cooled by collisions wi th 



protons having Tp > (lEsser et al. 



1999 



Landi &: Cranmer 



2009 



Hahn &: Savin 



20131 ) 



Thus, throughout the observed height range we expect Tj > Tg. We have applied this 
constraint to our analysis by rejecting fits th at imply Tj < 8 x 10^ K, which is consistent 



with a previo us analysis of this o bservation (IHahn et al. 



coronal hole (IWilhelm et al 



201ll ). 



2OI2I ) and is a typical Te for a 



4. Wave Amplitude and Energy 

Fits to equation ([8]) were performed over the range 1.02 - 1.12 Rq. The fits used a total 
of eleven lines from five ions (see column three of Table [1]). Figure |2] shows an example 
of the fit for Fe XI. If all the lines originate from the same volume, then it is expected 
that they will all have the same fnt, though not necessarily the same fth- As expected, 
the inferred fnt from the five different ions were in reasonable agreement with one another 
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(Table |2]). Thus, took the unweighted mean of the results from the different ions to find 
that v^t = 33.0 ± 2.4 kms'^ at 1.05 Rq. 



In the study of 



Dolla fc SolomonI (120081 ) they found v^t = 15 ± 2 kms ^, which is 



significantly smaller t han we find here. 



apparent discrepancy. 



here are several possible explanations for this 



Dolla fc SolomonI (120081 ) focussed on an observation made in May 



2002. This time period was near solar maximum, whereas our data were obtained near solar 
minimum and so the difference may reflect some solar cycle variation. Additionally, they 
described the polar coronal hole as "not well developed" and so their data likely contains 
other structures along the line of sight. For their analysis they studied a line from Mg x, 
which is a lithium-like ion with a peak formation Tg ~ 1.3 x 10^ K, but also h as a tail 



20091). Thus, if 



of high ion abundance towards much greater temperatures (iBryans et al. 
fnt is smaller in hotter structures than it is in a coronal hole, then their value would be 
systematically underestimated. 

Figure |3] shows function of R for Si vil 275.37 A, Fe ix 197.86 A, Fe x 184.54 A, 

and Fe xi 188.22 A. These lines, listed in the fourth column of Tabled] were chosen because 
they could be observed to relatively large heights. In each case the corresponding fth from 
Table |2] has been subtracted from Ves using the values for each ion determined from the fits 
(see Table E]). The solid line in the figure shows the mean fnt(-R) in 0.03 Rq bins for these 
lines and the dashed line shows the predicted ne trend for undamped waves. The data 
show that is consistent with undamped waves below about 1.15 Rq. However, at larger 
heights fnt deviates from this profile, which implies wave damping. 



The energy den s ity flux carri ed by the waves can be estimated using ( jPoyle et al. 



Banerjee et al. 



1998 



Moran 



20011 ) 



19981 : 



F = 2pvlVA, 
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where p is the mass density and 



is the Alfven speed with B being the magnetic field strength. To estimate the varying 
magnetic field strength for the supe rradially expan d ing po lar coronal hole we used the 



empirical model from equation (6) of 
A{R)/A{Rq) this gives 



Cranmer et al. 



( Il999bl ) . In terms of the area expansion 



B{R) = B{Rq)A{Rq)/A{R). 



Wangl (120 lOf ) gives B{Rq) = 7.3 G for this solar minimum. Other me asurements s . 



low that 



20091 ) ■ so 



the polar magnetic field can vary by a few Gauss between solar cycles ( iWang et al 
we estimate the uncertanity within a single solar minimum to be about 1 G. At low heights 
p can be found from the measured Ue- For larger heights it was necessary to extrapolate 



the de nsity measurements. We did this using the profile from 



Cranmer fc van Ballegooijen 



( I2OO5I ) which was based on white light measurements out to several Rq. Their ne{R) 
function was scaled to match our measurements at 1.12 Rq. The uncertainty of the scaling 
factor was taken to match that of rie at 1.12 Rq. Additionally, we found that F remains 
nearly identical if we simply use the hydrostatic fit for over the entire height range. 

Figure H] shows the energy density fiux F as a function of height based on the averaged 
results for fnt, plotted in Figure El These data are also listed in Table [3l The error bars 
represent the combined uncertainties from v^t, p and B. One can see that F is decreasing 
with height, but some of this decrease is due simply to the expansion of the coronal 
hole. The dashed line in Figure H] illustrates the variation of F with height for undamped 
waves, where the decre ase is due only to the superradial expansion of the coronal hole 



( ICranmer et al. 



1999b[ ). Clearly the waves are damped more rapidly with height than 



predicted by expansion alone. 



To more clearly show the effect of damping we show the quantity FA{R)/A{Rq) in 
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Figure [5] (also listed in Table [3]). In this plot, measurements for undamped waves would 
fall on a horizontal line. The data are consistent with undamped waves for i? < 1.15 -R0. 
The dashed line in Figure is drawn at the average of the points below 1.12 Rq, which 
is F = 6.7 ± 0.7 X 10^ erg cm~^s~^. This is the ani ount of Alfven wave energy present 



at the base of the corona. 



Withbroe Sz Noyed (jl977l ) estimated that 8x10^ erg cm ^ s ^ 



is typically required to heat a coronal hole and accelerate the fast solar wind. However, 



during the 2007 - 2009 solar minimum the solar wind was observed to be unusua 



being about 25% less energetic than in previous minima (iMcComas et al. 



2008 



ly wea k, 



Wang 



20101 ). This implies that for the recent solar minima only about 6 x 10^ ergcm"^ s~^ would 
be required. Our measurements show not just that the amount of energy carried by the 
waves is sufficient to account for coronal heating and solar wind acceleration within the 
coronal hole. In addition they indicate that the waves are indeed damped, losing ~ 85% of 
their initial energy 1.5 Rq (Figure [5]). Thus, our findings indicate that the waves do in fact 
provide most of the required heating. 

The length and time scales over which the waves are damped provide benchmarks for 
theoretical calculations. In order to estimate the length scale over which the waves are 
damped we fit an exponential to FA{R)/A{Rq). This fit is illustrated by the solid line in 
Figure The initial value of F, the height where damping begins Rd, and the exponential 
damping length L^, were free parameters of the fit. The initial F was the same as found 
above, F = 6.7 ± 0.7 x 10^ ergcm'^s'^ The fit found Rd = 1.12 ± 0.04 Rq, which is 
consistent with the point where v^t deviates from the Ue^^'^ trend (e.g.. Figure [3]). The 
relatively large error bar is due to the course binning used here. We find that the damping 
length is about Ld = 0.18 ± 0.04 Rq. We can also estimate a timescale for the damping. 
Since the velocity of the waves is about the Alfven speed, we can convert R to wave travel 
time t. Fitting an exponential to the data as a function of t, we find that the damping time 
is about 68 ± 15 s. 
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Solar wind models show that in order to accelerate the fast solar wind to the speeds 
observed far from the Sun, some in put of wave ene rgy is needed above the point where 



the solar wind becomes supersonic (jCranmer 



2002I ). This suggests that not all of the 



wave energy should be damped at very low heights. The amount of initial energy that 
is required to be undamped to large heights is about 1 x 10^ ergcm~^s^^. Our results 
show that at least up to about 1.4 Rq sufficient energy remains in the waves to provide 
the additional acceleration for the solar wind. However, because of the large uncertainties, 
our data can also be consistent with the wave energy going to zero at large distances. 
To get a rough estimate we performed a similar fit to the one described in the above 
paragraph, but using an exponential plus a constant. We find that at large distances 
FA{K)/A{Rq) 0.6 ± 1.4 X 10^ ergcm~^s~^. One additional source of uncertainty for 
this estimate is that we do not know how the ion temperatures are changing with height. 
As the dissipation of the wave energy is likely to heat the ions, the assumption of constant 
temperature is probably becomes less reasonable at the larger heights in our observation. 
Since any resulting thermal broadening would increase with height, our assumption of 
constant Vth would cause us to underestimate the actual damping. 



5. Ion Temperatures 

The temperature of each ion can be inferred at 1.05 Rq from all the line widths 
observed at that height (column five of Table [1]) by subtracting the non-thermal width 
fnt = 33.0 ± 2.4 kms~^. The circles in Figure [6] shows Ti for each of the ions measured as a 
function of charge to mass ratio q/M, in units of elementary charge e per atomic mass unit 
(amu). These data are also given in Table |H 

Previous measurements have found that Ti is greater than To for q/ M < 0.2, while for 



slightly higher q/M ions Ti ^ T^, but may increase again for q/M > 0.3 ( iLandi fc Cranmer 
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20091 : iHahn et al.ll2010l ). Here we find a similar pattern with respect to q/M. For q/M < 0.2 
Tj ?5i 2 X 10^ K, while for larger q/M T; 1 x 10^ K, which is about the expected value of 
Tg for a coronal hole. For larger q/M > 0.3 there is a suggestion that increases if one 
looks at the S x and O vi data. The yet higher q/M point from Si x might appear to 
contradict this trend, but this could be a systematic error since Si x is formed at relatively 
high temperature s and so a large fr action of the emission may come from structures outside 



the coronal hole (IHahn et al 



2OI2I ). We should also note that our uncertainties are large 



enough that we cannot rule out that Ti is actually constant over the entire range with 
respect to q/M. 

Different models for ion heating predict different dependence s of Tj on q/M . Thus, 



Ti measurements can be used to test these models. For example. 



Cranmer et al. 



fll999ah 



developed a model in which the ions are heated by ion cyclotron waves. Based on their 



equations (2) and (15^ 



spectral index of 3/2 (ILeamon et al. 



, ignoring col 



isions , and assuming a typical solar wind p lasma wave 



1998 



Podesta et al 



2007 



Chandran 



2010) one finds 



r,ocM|-l)"%i 



M 



(9) 



Another model is that of 



ChandranI (120101 ) in which Alfven wave turbulence causes ion 
orbits to become stochastic and absorb energy from the turbulence. They derive a 
dependence of T[ on q/M if (a) the ratio of the turbulent velocity fluctuations to the thermal 
velocity perpendicular to the magnetic field are the same for all the ions, (b) the turbulent 
fluctuations have spectral index 3/2, and (c) collisions can be ignored, finding that 



(10) 



However, the neglect of collisions in deriving both of these trends is probably not a very 
good approximation at these low heights. Nevertheless, we can compare these predictions 
to our data. The open squares and diamonds in Figure [6] illustrate the predicted trends 
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from 



Cranmer et al 



fll999ah and 



Chandraru (|2010| ). respectively. In each case the theoretical 



trends have been multiplied by a scaling factor that was chosen to produce the best average 
agreement with the observations. Given the large uncertainties in our analysis and the 
neglect of collisions in the models, both predictions show reasonable agreement with our 
data. 



6. Summary 

We have shown that Alfven waves in a polar coronal hole possess sufficient energy 
to heat the coronal hole and that this energy is actually dissipated from the waves 
at sufficiently low heights. To show this we determined separately the thermal and 
non-thermal components of spectral line broadening in a coronal hole. Our method relies 
on the observation that waves are undamped at very low heights and on the assumption 
that the temperature of each ion does not change with height at low heights. From the 
derived v^t "we show that the energy carried by the waves is 6.7 ± 0.7 x 10^ ergcm~^s^^, 
which is sufficient to heat the coronal hole and accelerate the fast solar wind. The waves 
are damped over a length of about 0.18 ± 0.04 Rq, with a corresponding timescale of about 
68 ± 15 s. Although our measurements are limited to i? < 1.5 Rq, they suggest that enough 
energy remains in the waves to provide the extended heating of the solar wind above the 
sonic point that models show is required. Additionally, we measured T[ for each ion to be 
in the range of about 1-2 MK. We found a weak trend where low q/M < 0.2 ions have 
the highest temperature, q/M ^ 0.2 - 0.25 are lower with T; ^ Tc, and q/M > 0.25 have a 
slightly increasing temperature. Our uncertainties are not sufficiently small to distinguish 
between the predictions of two ion heating models. Those models are also not realistic for 
these heights since they neglect Coulomb collisions. Our results, though, do demonstrate 
that such a comparison is possible in principle, with some additions to the model and higher 
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Fig. 1. — The filled circles indicate the electron density derived from an Fe ix inten- 
sity ratio. The solid line shows the hydrostatic equilibrium fit to the data using equa- 
tion (171) in the range 1.02 - L12 i?©. The dotted line shows the empirical model from 
Cranmer &: van Ballegooijenl ( 120051 ). scaled to match the data at 1.12 Rq. For the analy- 
sis we used the hydrostatic fit for R < 1.12 Rq and extended to larger heights using the 
empirical model. 
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Fig. 2. — The filled circles show the measured effective velocity Ves for Fe xi 188.22 A. The 
solid line illustrates the average fit to the data between 1.02 and 1.12 Rq using equation ([8]). 
The fit parameters for this ion were Vth = 25.8 ± 5.4 kms~^ and fnt = 32.2 ± 4.2 kms~^. 
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Fig. 3. — Symbols indicate the non-thermal velocity v^t from the strongest observed lines. 
The filled circles and solid line show the averaged v^t combining the data from the various 
ions. The dashed line illustrates the predicted electron density rie^^^ trend for undamped 
waves. 
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Distance from Sun Center (Rq) 



Fig. 4. — Wave energy density flux F as a function of height (filled circles). The dashed line 
illustrates the predicted trend for undamped waves. 
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Fig. 5. — Points indicate F multiplied by the expansion factor A[R)/A{Rq) in order to more 
clearly show the decrease due to damping. For undamped waves, the points would fall on 
a horizontal hne. The dashed line indicates the average of the points below 1.12 Rq. This 
average shows that F = 6.7 ± 0.7 x 10^ ergcm~^s~^ is present in the waves at 1 Rq. The 
solid line gives an exponential fit, from which a damping length of 0.18 Rq was derived. The 
point at which the exponential decay begins was a free parameter of the fit, with the result 
= 1.12 ±0.04 
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Fig. 6. — Filled circles show the temperature for each ion, derived by subtracting the 
average fnt = 33.0 ± 2.4 kms~^ from Ves, plotted versus charge to mass ratio q/M for 
different ion species. The open squares and open diam onds show the patter n of Ti versus 



q/M predicted by the models of 



Cranmer et al. 



(Il999al ) and IChandranI (l201Cll ). respectively. 
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Table 1:: Line List. 



Ion 


A [A) 


tjq. 


Used for: 
til v^t{n) iil^i.Uo 


U IV 


z ^ y.ooi 




* 


U IV 






* 


Mg VII 


Z /o.i54 




* 


bl VII 


O TO /I Q 

z / z.o4o 


* 


* 


bl VII 


Z / b.ooi 


* 


* * 


bl VII 


Z ^ 0.0 


* 


* 


bl X 


Zoo. 6 (4: 




* 


bl X 


zoi.Uo/^ 




* 


bl X 


z i 1 .yyz 




* 


bl X 


z / / .zo4 




* 


Fe VIII 


ioo.zio 


* 


* 


Fe VIII 


iob.oyy 


* 


* 


Fe VIII 


194.661 


* 


* 


Fe IX 


188.497 


* 


* 


Fe IX 


1 on 0/11 


* 


* 


Fe IX 


197.862 


* 


* * 


Fe X 


184.537 


* 


* * 


Fe X 


190.037 




* 


Fe X 


193.715 
257.259 




* 
* 


Fe X 


257.263 




* 


Fe XI 


180.401 




* 
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Ion A (A)^ Used for: 

Eq. iD Fit Vnt{R) Ti(l.O5i?0) 
Fe XI 188.217 * * * 

Fe XI 188.299 • • • ^ 



^Wavelengths from CHIANTI flDere et al. 



1997 



Landi et al. 



20121). 



^AA was constrained to be identical for Fe xi 188.217 A and 188.299 A. 
^Brackets indicate blends from the same ion. 
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Table 2:: Values for t>th and fnt at 1.05 Rq from fitting 
Equation ([8]) over 1.02 - 1.12 Rq. 



Ion 


Vth (kms"^) 


Vnt (kms"^) 


Si VII 


23.5 ± 1.5 


33.6 ± 1.2 


Fe VIII 


19.9 ±2.7 


29.8 ± 1.8 


Fe IX 


20.4 ±3.0 


34.9 ±1.8 


Fe X 


18.7 ±2.6 


34.5 ±1.5 


Fe XI 


25.8 ±5.4 


32.2 ±4.2 
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Table 3:: Non-thermal velocity and Energy Flux Density. 





^T' ^ 1 n^ OTfT r*m ^ c ^\ 

r I lu t;!^ Cili o J 




eig cm b j 


1 no 


0.0 IC l.O 


6.1 


±1.4 


i.UO 


O.D It 1.1 


7.1 


±1.4 


l.Uo 


A -I- n Q1 
4.oo zt u.yi 


6.5 


±1.3 


111 
i.ii 


/I /I K -L n Q"? 
4.40 zt u.yo 


7.3 


±1.5 


11/1 

1.14 


Q 71^ -I- n 
o. / zt u.yo 


6.9 


±1.7 


117 
1 . 1 ( 


Qn -1- n c;7 

Z.OU It U.O ( 


4.7 


± 1.2 


1 on 

l.ZU 


o 1 /I -I- n sc; 

Z.14 Zt U.oO 


4.9 


± 1.9 


1 0"? 
1 .zo 


1 o"^ -1- n 00 
i.zo zt u.zy 


3.18 ±0.73 


1.26 


0.91 ±0.39 


2.6 


± 1.1 


1.29 


0.89 ±0.33 


2.8 


±1.0 


1.32 


0.64 ±0.16 


2.20 ±0.57 


1.35 


0.65 ± 0.42 


2.5 


±1.6 


1.38 


0.48 ±0.16 


1.98 


±0.68 


1.41 


0.34 ± 0.33 


1.5 


±1.5 


1.44 


0.20 ±0.10 


0.98 ±0.50 
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Table 4:: Ion Temperatures. 



Ion 


q/M (^) 

^ ' ^ amu ' 


Ti (MK) 


VI 


0.31 


2.9 ±1.2 


Mg VII 


0.25 


0.91 ±0.44 


Si VII 


0.21 


1.04 ±0.47 


Si X 


0.32 


1.41 ±0.45 


S X 


0.28 


1.44 ± 0.57 


Fe VIII 


0.13 


1.91 ±0.74 


Fe IX 


0.14 


2.32 ±0.86 


Fe X 


0.16 


1.50 ±0.96 


Fe XI 


0.18 


2.74 ± 0.94 
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